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Second and first order phase transition in 
three dimension Gross-Neveu mo de£| 

Zhou Bang-Rong 

Department of Physics, Graduate School of the Chinese Academy of Sciences, Beijing 100039 

Symmetry restoring phase transitions in three dimension Gross-Neveu model are shown to be 
second order at finite temperature T and first order at T = and finite chemical potential jj, by 
critical analysis of the dynamical fermion mass based on the gap equation. The latter is further 
verified by effective potential analysis. The resulting tricritical point is (T, /i) = (0,m(0)), where 
m(0) is the dynamical fermion mass at T = /x = 0. Physical difference between the above second 
and first order phase transition is illustrated by means of variations of thermodynamical particle 
density. 

Keywords: thermal gap equation and effective potential, symmetry breaking and restoration, second and 
first order phase transition, tricritical point, particle density 



The Gross-Neveu (GN) model Q with four- 
fermion interactions is a good laboratory to re- 
search symmetry restoring phase transitions at 
high temperature and high density H, IE ^ Hi- 
One believes that some phase structure of Quan- 
tum Chromodynamics (QCD), e.g. a tricriti- 
cal point in temperature T-chemical potential fj, 
phase diagram, i.e. the point where second order 
and first order phase transition meets, could ap- 
pear in such type of models Existence of a 
tricritical point is of both theoretical and exper- 
imental significance. Since the three dimension 
(3D) GN model is easy to deal with analytically, 
and in addition, in this model there is still a ques- 
tion to be clarified further of that whether first 
order phase transition exists or not in it |^ , 
it is interesting to make a careful analysis of the 
phase transition of this model at finite tempera- 
ture, especially at high density. In this paper, 
we will report our results and indicate that a 
tricritical point exists indeed in 3D GN model. 
Our method is to analyze critical behaviors of the 
dynamical fermion mass as the order parameter 
of symmetry breaking based on the gap equa- 
tion obeyed by the dynamical fermion mass and 
on the corresponding effective potential at finite 
temperature and chemical potential. The discus- 
sions will be made in the fermion bubble diagram 
approximation which is equivalent to the leading 
order of 1/N expansion. 

We first briefly review the used model l7l| . The 
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Lagrangian will be written by 
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k=l 



(1) 

where ^pk{x) is a spinor with N "color" compo- 
nents. In 3 dimensions, the coupling constant g 
has mass dimension -1 thus the model is pertur- 
batively non-renormalizable however it is renor- 
malizable in 1/N expansion Q. We take the 
7^^ [jj, = 0,1,2) to be 2 X 2 matrices and re- 
late them to the Pauli matrices cr*(i = 1,2,3) 
by 7° = CT'^, 7^ = ia^ and 7^ = ia^ . Noting 
that in three dimensions there is no "75" matrix 
which anticommutes with all the three 7^ matri- 
ces. Consequently C{x) has no chiral symmetry. 
However, the action of the system is invariant un- 
der the special Vi and 1^2 parity transformation 
(omitting the index k of ip) 



tp{t,x\x^) ^ j^i;{t,x\-x^) 



and the time reversal 



■>p(t,x\x'^) 7^V(-^,a;\x2). 

If the four- fermion interactions in Eq.(l) will 
lead to the dynamical mass term —rmpip, then 
the discrete symmetries Vi, V2 and T will be 
spontaneously broken. These broken symme- 
tries, however, could be restored at high tem- 
perature and/or high density. This is just the 
problem we will discussed. 

At zero-temperature T — and zero chemical 
potential /i — 0, assume the four-fermion interac- 
tions gJ2^=i (''/''^V'fc) /2 can lead to the fermion 
condensates 0; then we will ob- 

tain the gap equation obeyed by the dynamical 
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fermion mass m(0) 



where 



(27r)^ P - m2(0) + ie 



(2) i^2(r,Af,TO) = Tin 



After the Wick rotation, angular integration and 
introduction of 3-D Euclidean momentum cut-off 
A, Eq.(2) becomes 



1 = 



gNK 



, m(0) A 
1 : — arctan 



A 



m(0) 



(3) 



For sufficiently strong coupling gN that 
gNA/n"^ > 1, Eq.(3) can be satisfied. When 
T 7^ 0, we must replace the vacuum expectation 
value X]fcLi(V''^'0fe) by the thermal expectation 
value '}Zk^i{w'4'k)T- In the real-time formal- 
ism of thermal field theory, this means the 
substitution of the fermion propagator Q 



- 27r(5(/2 - m2) 
(4) 



/ — m(0) -I- ie 
with 

sin2 6'(/°,/i) 



I — m + ie 



9{P) 



exp[/3(?o - yu)] + 1 



(5) 

where (3 = 1/T and m = m(T,ii) is the dynam- 
ical fermion mass at finite temperature T and 
finite chemical potential /z. As a result, the gap 
equation at T 7^ becomes 



gNA 



m A TT 

1 - — arctan —F2 (T, /x, m 

A m 2A 



(6) 



1 + e-(™-^)/^J + i-fi ^ fi). 

(7) 

Considering the gap equation (3) at T = and 
the fact that the momentum cut-off A may be 
high enough so that arctan 

we will reduce Eq.(6) to 



A 

m(0) 



arctan ■ 



to(0) = m + F2 (T, n^m). 



(8) 



It is easy to verify that dF2{T, fi, m)/dT > and 
dF2{T, ^,m)/dfi > 0, i.e. F2{T,fi,m) is respec- 
tively a monotone increasing function of temper- 
ature T and chemical potential /i. As a result, 
it can be seen from Eq.(8) that, as T and/or /x 
increase, m must decrease and finally go to zero 
at a critical point (TcHc)- The critical equation 
to determine the point {Tc,^c) will be 



m(0) = F2(Te,/ic,m = 0) 



(9) 



In the above discussions, we have actually used 
an implicit assumption, i.e. m could go to zero 
continuously as T and/or /i increase. This just 
represents the characteristic of a second order 
phase transition. Hence the T — /i phase diagram 
given by Eq.(9) will be the one of second order 
phase transition. The critical behavior of m near 
a critical point can be derived from Eqs.(8) and 
(9). Since near a critical point (Tcfic), m « 0, 
thus m/T <C 1 for a finite temperature T. Then 
Eq.(8) with Eq. (7) will lead to 



i(0) = Tin (l + e'"/^) + Tin (^1 + e"^/^) + 



2r[l + cosh(/x/r)] 



o 



m 
T3 



(10) 



Substituting Eq.(9) with /Xc replaced by /i into Eq.(lO) we will obtain the critical behavior of near 
Tc when /i is fixed 



2r(l-Hcosh|;) l\n 



2(1 + cosh ^ 

J r. 



|-tanh^|(T,-r) whenr<r, (11 



Substituting Eq.(9) with Tc replaced by T into 
Eq.(lO) we will obtain the critical behavior of 
m? near /x^ when T is fixed 

m2 = 2Tsinh(/VT)(^c - A^), 

when ^J.< ficandT ^ 0. (12) 



Eqs. (11) and (12) show that when T ^ 0, the 
symmetry restoring phase transitions whether at 
Tj, or at fic are second order. However, the above 
analyses are inapplicable to the case of T = 0, 
since the approximate expansion (10) is obviously 
not valid in this case. 
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For T — 0, we must come back to the gap 
equation (8). Considering that 

, N / when ^ < m 

lim i'2{i , U, m) = < , 
T^o ^ ' ' I /i — TO when fj, > m 

(13) 

we can obtain from the T ^ hmit of Eq.(8) 
that 

to(0) — TO, when ii < m (14) 
to(0) — /i, when /i > m (15) 

Eq.(14) can be changed into that m = to(0), 
when fj, < m = to(0). Then Eq.(15) wiU become 
to(0) = /i, when > to(0), two contradictory 
formulas each other. This implies that the gap 
equation can not be satisfied when > to(0) thus 
there is no symmetry breaking and the order pa- 
rameter m must be zero. The above result can 
be summarized as 

^ ^ f to(0) when < m(0) 

1 when /i > to(0) ' ^ ' 

Eq.(16) indicates that when T ^ 0, the order 
parameter to will jump from m(0) down to 
as /i crosses over the critical chemical potential 
/Zc = to(0), hence the phase transition is first or- 
der. This conclusion can also be verified by an 
effective potential analysis. It can be proven that 
the gap equation (8) may also come from the 
extreme value condition of the effective poten- 
tial v}fj{T,fi,m), i.e. dv}f'>fiT,fi,m)/dm = 0, 
where 

^"g^ = £ _ ^(0) + F,{T, m)] , 

(17) 

Thus we can obtain from Eq.(17) by integration 

V}fl{T,fi,m) = — dm'm' [to' - m(0) 
27r Jo 



+F2{T,fi,m')] 



(18) 



which has been normalized to V^^ ^{T, ii,m — 
0) = 0. By means of Eq.(13) we may obtain 
the T hmit of V^f^{T, fi, to) 



V^ffiT = 0, M, m) = — {3to2[^0(^ - to) 
-m(0)] -I- e{m - fi){2m^ + f^^)} ■ 



(19) 



The shape of V^^j{T = 0,/i, m) will vary as fi 
increases. 



1) /i = 0. V^^j|(T — 0, fi, m) will have the maxi- 
mum point TO = and the minimum point to = 
to(0) with v}fj[T = 0,M,m(0)] = -^^^^^(O). 
This simply represents the spontaneous symme- 
try breaking at T = /i = 0. 

2) ^ < to(0). V^^j{T — 0,/i, to) has the same 
maximum point to = and the same minimum 
point TO = to(0), but the minimum V^j:j[T = 

0,/x, to = to(0)] = jj^ifJ-^ — TO'^(O)] will rise as n 
increases. 

3) n = to(0). The effective potential now be- 
comes 



K5}[T-0,Ai = TO(0),TO] 

= -^0[m ~ to(0)] [to^(O) - 3to2to(0) + 2to^] 



which indicates that the total real axis segment 
< TO < to(0) including the original to = to(0) 
when /i < to(0) will be the minimum points of 
V^f\[T^O,^l = m{0),m]. 



4) fi> to(0). v}fj[T = 0, /i > to(0), to] will only 
have a minimum point to = and this implies 
that the broken symmetries will be restored. Ob- 
viously, the critical chemical potential iic should 
be determined by the condition 



0,/i, to) 



dm? 



to(0) = 0, 



m=0 



i.e. at /ic, TO = must change from a maximum 
point into a minimum point. Eq.(20) gives that 
/Xc = to(0). As /X increases and crosses over to(0), 
the minimum point to of V^j^{T = 0,/i, to) will 
jump from to(0) down to 0. This will reproduce 
Eq.(16) and verify that (T, /i) = (0,to(0)) is a 
first order phase transition point. On the other 
hand, it is noted that (T, /i) — (0, to(0)) is also a 
solution of the critical equation (9) correspond- 
ing to second order phase transition when we take 
fic = 1^ and the limit Tc = T ^ 0. Hence we can 
reasonably conclude that (T, /i) — (0,to(0)) will 
be a tricritical point in the T — /i phase diagram 
of 3D GN model. 

The physical difference between the above sec- 
ond order and first order phase transition can 
be clearly displayed by variations of thermody- 
namical particle density near critical points. We 
will respectively consider the phase transitions at 
T ^ and T = 0. The general total density of 
fermions with mass to in 3 dimensions may be 
expressed by 



4 



f P 



m 
2^ 



dz z 



1 

1 



gy^-r + 1 



(— r — > r) 



, y = /3m, r = 



Tm , 1 + 6-!'+'' 

In 

2tt 1 + e-2'-'' 



— V ^-4 sinh(A;r)e-'=^ 

fe=i 



(21) 



For the model of dynamical spontaneous symmetry breaking, m should be identified with the dynami- 
cal fermion mass. When T 7^ 0, the phase transition is second order thus near a critical point (Tc, /Xc), 
y 1, we can make the expansion of n{T, /i, m) in the powers of y up to 0{y^) order and obtain 



n{T, fi, m) ~ n(T, /x, m = 0) 



e'' 



1 , m , 

-, when — <gC 1. 

47r e'- + 1 ' T 



(22) 



Now if we consider the phase transition at finite chemical potential and T ^ 0, then can get by using 
Eq.(12) 



n(T, /U, m) 



n{T, ii,m = 0)- 2}{e-+i) (si'^^ t) (/^c ~ /^)' "^^^'^ M ~ /"c 
n(T, n,m = 0), when n > Hc 



However, 



dn{T, jj, m) 
dfj, 



+ 27r(e'- + l) ^^^'^ T' WneU ^ ^ic 



dn{T,fj.,m=0) 

V 



when ii> He 



(23) 



(24) 



Eqs.(23) and (24) indicate that when /i ^ 
/Uc, n{T,iJ,m) will increase continuously up to Eq.(25) and obtain 
n{T,H,m = 0) but dn{T,ij,,m)/d^ will have a 
discontinuous jump across over /i^- These just 

show second order feature of the phase transition. n{T — 0, /j., m) = 

For discussion of the phase transition at T = 0, 
we first find out the T — > limit of n{T,n,m) 
from Eq.(21) 



((0) is first order, we can substitute Eq.(16) into 



0, when /x < m(0) 

47r' 



^, when /X > m(0) 



(26) 



n(T = 0, /X, m) = ^(/x — m) 



/x^ - 



(25) 



Eq.(26) shows that when T = 0, as /x crosses 
over the critical point /Xc = to(0), the particle 
density will jump from to /x^/47r (the density 
of massless fermions) and this just represents a 
characteristic of first order phase transition. 



Since the phase transition at T = and fi — He 
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